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Abstract. The generation of laser-driven dense relativistic electron layers from ultra-thin foils and their use 
for coherent Thomson backscattering is discussed, applying analytic theory and one-dimensional particle- 
in-cell simulation. The blow-out regime is explored in which all foil electrons are separated from ions by 
direct laser action. The electrons follow the light wave close to its leading front. Single electron solutions 
are applied to initial acceleration, phase switching, and second-stage boosting. Coherently reflected light 
shows Doppler-shifted spectra, chirped over several octaves. The Doppler shift is found oc 7^; = 1/(1 — 0$.), 
where f3 x is the electron velocity component in normal direction of the electron layer which is also the 
direction of the driving laser pulse. Due to transverse electron momentum p y , the Doppler shift by 4-y% = 
4j 2 /(l + (py/mc) 2 ) « 27 is significantly smaller than full shift of Aj 2 . Methods to turn p y — > and to 
recover the full Doppler shift are proposed and verified by ID-PIC simulation. These methods open new 
ways to design intense single attosecond pulses. 

PACS. 41.75.Jv Laser-driven acceleration, - 52.38.-f Intense particle beams and radiation sources in 
physics of plasmas, - 52.59.Ye Plasma devices for generation of coherent radiation 



1 Introduction 



There is a quest for high-quality x-ray sources in many 
fields of science [TJ. Coherent VUV- and X-ray sources 
now become available at XFEL facilities \2\ ■ Also the gen- 
eration of high harmonics from optical laser pulses has 
provided us with very useful VUV-sources recently [3j- 

In this paper, we deal with Thomson scattering from 
relativistic electrons. Incoherent Thomson scattering con- 
verts photons of incident frequency wj, to frequencies lo = 
-f 2 (l — (3cos9) 2 u)L due to relativistic Doppler shift [4]; here 
[3 = v/c is electron velocity normalized to light velocity 

c, 9 is the scattering angle, and 7 = l/yl — (3 2 the usual 
relativistic factor. Thomson scattering has been measured 
from high-7 electrons at conventional accelerator facilities 
jS] and recently also from laser generated electron beams 
[6]. In these cases, light was backscattered incoherently. 

Coherently backscattered photons can be obtained from 
electron layers provided they are dense enough so that 
several backscattering electrons reside in the volume A^j , 
where A r is the wavelength of the light in the rest frame 
of the electron layer. When moving close to velocity of 
light, we refer to these layers as relativistic mirrors with 
7-factors 7 m [7]. Such mirrors have been created recently 
in form of nonlinear plasma waves near the point of wave- 
breaking when they develop diverging density spikes; val- 
ues of 7,„ =4—5 could be deduced from the backscattered 
spectrum [8]. 



Relativistic mirrors are also formed when irradiating 
solid surfaces with laser amplitudes ao = eA^jrac 2 > 1 
[9|10|llj : here e and m are electron charge and mass, re- 
spectively, and Alo is the amplitude of the vector poten- 
tial of the light wave. For ao > 1, the solid surface turns 
into dense plasma, and electrons move close to velocity of 
light. For linear polarized light, the ponderomotive force 
has a component oscillating at twice the laser frequency 
wj,. This leads to an oscillating mirror front and to high 
harmonics in the reflected light. These surface harmonics 
have been observed recently [12] • The results show har- 
monic orders up to 3200 corresponding to 3.7 keV pho- 
tons. They are in agreement with theory [13] . 

In the present paper, we consider laser irradiation of 
ultra-thin (nm) foils [14115116117) at intensities high enough 
to separate all electrons from ions. The electrons then 
move with the front of the laser pulse, forming a dense 
relativistic electron layer. This regime has been investi- 
gated recently by Kulagin et al. [18J . Here it is studied for 
Thomson backscattering of a counter-propagating probe 
pulse. 

For separating electrons from ions, the laser field ao = 
Elq has to be larger than the total electric field £0 that 
builds up due to charge separation. Here both the laser 
field Elq = k^A^o with = lul/c and the longitudinal 
field £0 ar e taken in units of Eq — mcLOi^je. For a pla- 
nar foil with electron density n e and thickness d, we find 
£0 = Nktd. Here N — n e /n crit = w 2 /w 2 L is the den- 
sity normalized to the critical density n cr i t and can be 



2 



J. Meyer-ter-Vehn and H.-C. Wu: Coherent Thomson backscattering from relativistic ultra-thin electron layers 



20 
15 
10 
5 

0, 



t/z L = 



t/T L =l 



t/T L =2 






x/X 



2 -2 



L 





xA,, 



2 -2 





xA, T 



Fig. 1. Irradiation of a graphene layer with a two-cycle laser pulse of ao — 5 (2D-PIC simulation). Interfaces of electron layers 
are plotted in the plane of laser polarization before as well as one and two laser cycles after the light front has touched the foil. 
The simulated laser intensity is shown in shaded grey; the broken line shows the immobile ion layer left behind. 



conveniently expressed by the ratio of plasma frequency 
uj p = (47re 2 n e /m) 1 / 2 and laser frequency lol. The condi- 
tion for complete electron blow-out is 

a Q > Nk L d. (1) 

As a reference foil material, we consider in the following 
graphene which consists of a mono-atomic layer of carbon 
atoms arranged in regular hexagons with a side length 
of s = 0.141 nm. Having an area of A — (3\/3/2)s 2 = 

5.2 x 10~ 16 cm 2 and a total of 12 electrons attached to 
each hexagon, the electron areal density is n e d = 12 /A = 

2.3 x 10 16 cm -2 , and one finds £o = Nk^d ~ 1. Graphene 
layers have been successfully produced in pieces up to 
100 x 100 fim 2 [T^] • This is sufficient in principle for experi- 
ments with existing laser pulses in the range of a w 3 — 10, 
corresponding to intensities of 10 19 — 10 20 W/cm 2 . It is 
interesting to notice that an area of 10 x 10 /itm 2 contains a 
total charge of about 1 nC, comparable and exceeding val- 
ues in other recent experiments using laser wakefield accel- 
eration and gas targets. Alternatively, multi-layer carbon 
foils with thicknesses of a few nm are available for exper- 
iments, but require higher laser intensities for complete 
electron blow-out. 

For illustration, a two-dimensional (2D) particle-in- 
cell (PIC) simulation of a graphene layer irradiated by 
a laser pulse is shown in Fig. 1. Electron density is plot- 
ted for three different times, and also the driving laser 
pulse is shown. Here we have chosen a two-cycle pulse 
E(t)/Eq = aosinr with < r < 47r and ao = 5. At 
the corresponding intensity of 3 x 10 19 W/cm 2 , all elec- 
trons within the focal spot are taken as free initially, i.e. 
ionization is not treated explicitly. It is seen that the elec- 
trons of the irradiated spot have separated from the ions 
and are moving with the laser pulse, surfing as a dense 
layer near the pulse front. One should notice that these 
ultra-thin foils are transparent to the light, even though 
they have densities far above the critical density. This 



is because they are much thinner than the skin depth 
A s fcm ~ c/ujp w Al/10. The light is therefore penetrating 
through the foil, and each electron is accelerated accord- 
ing to the local laser field. This means that we can use the 
analytic solution for relativistic electron motion in a plane 
laser field as an approximation, supplemented by the space 
charge electric field E x . It is constant in time for each sub- 
layer of the foil, as long as these sub-layers keep their rela- 
tive order and do not overtake each other. This option for 
analytical treatment is an important feature of ultra-thin 
foil acceleration. It is quite different from usual relativistic 
laser plasma interactions, e.g. encountered in laser wake 
field acceleration, where often complicated plasma fields 
dominate. 



Here it is direct interaction with the external laser 
fields that accelerates the electrons. As we shall discuss 
in section 2, electrons gain momentum up to a maximum 
of p m ax = {Px,P v ) = Qa\mc, 2a mc) and energy up to 
W max = (7 — l)mc 2 = 2a 2 >mc 2 in the first half-wave of 
the light pulse, at least those at the layer front where 
E x = 0. For the laser pulse ao — 5 used in Fig.l, the 
electron energy rises up to W max « 25 MeV according 
to these estimates, but then drops again to zero in the 
second half of the laser wave. This behavior reflects the 
Lawson- Woodward theorem [20], stating that there is no 
net energy transfer to charged particles by light waves in 
vacuum. Having in mind generation of X-rays by Thom- 
son scattering, this is acceptable, because there is plenty of 
time for back-scattering of a probe pulse by the dense elec- 
tron sheath while surfing on the wave. Actually, the time 
available is approximately t = (3ag/8)Ti which is about 
10 laser periods T£ for a — 5. For 7 « 2a 2 , = 50, one 
might be tempted to expect coherent photons Doppler- 
shifted by a factor 47 2 = 10 4 . Unfortunately, this is not 
the case because the effective 7-factor at which the layer 
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is moving in normal ir-direction is 



lx = = ll\J l + {Py/mcf « (2) 

and the frequency is shifted by a factor w 2j rather 
than 4j 2 . The reason for this is inherent to direct laser ac- 
celeration which necessarily involves large transverse mo- 
menta p y oc eto for large 7 oc p x oc ag. Electrons move 
under an angle « p y /p x oc 1/ao relative to the laser di- 
rection. For incoherent Thomson scattering, the reflected 
signal is emitted close to the direction of the individual 
high-7 electrons and show the angle 0. In the case of co- 
herent scattering from a planar layer, however, this main 
component is suppressed by destructive interference and 
only the component scattered into the normal direction of 
the sheet will survive. Since incoherent scattering is pro- 
portional to the number N s of scatterers, while coherent 
scattering is oc N%, the latter is bound to win as the layer 
becomes dense enough. In practice, one may encounter 
cases intermediate between incoherent and coherent scat- 
tering and observe the change in direction, frequency, and 
intensity. 

The present paper is intended to explore basic fea- 
tures of electron sheet acceleration and back-scattering in 
ID geometry. The analysis is restricted to schematic laser 
pulses of single-cycle and half-cycle form. In section 2, we 
review the analytic solution of single-electron motion in 
a planar light wave and discuss three relevant situations: 

(1) acceleration of an electron initially at rest, (2) collision 
of a relativistic electron with a counter-propagating wave, 
and (3) boost acceleration of an electron entering the light 
wave already at relativistic energy. Results of cases (1) and 

(2) are then compared with ID-PIC simulations in sec- 
tion 3, including back-scattering of a counter-propagating 
probe pulse. The coherent reflectivity of an ultra-thin rela- 
tivistic layer is derived analytically in a companion paper, 
referred to as paper II [22]. The full analytical descrip- 
tion of the driven electron layer, including self-radiation 
and space charge fields, was developed by Kulagin [23] : 
it will be applied to the present reference case in another 
companion paper [24]. 



2 Single electron in plane laser wave 

2.1 Equations of motion 

Here we review the analytical treatment of a relativistic 
electron interacting with a plane light wave in vacuum, 
to the extent it is needed in the following discussion [25] . 
The linear polarized light wave propagates in x-direction 
and has electric and magnetic field components E y and 
B z , respectively. In addition, we allow for a longitudinal 
E x field. The equation of motion then read 



dp x 
dt 
dp y 
~~dt 



-E x — (3 y B z , 
-E y + fl x B z , 



7 2 = 1 +Pl +P 2 y , 

-rr = —PxE x — (3 y E y , (3) 

-jr = Px =p x ft, 
at 

d y a 1 

Here we use dimensionless variables corresponding to di- 
mensional ones according to 

(t,x,y) = (uj L t,k L x,k L y), 
(Px,P y ) = (px/mc,p y /mc), 
(E x ,E y ,B z ) = (E x ,E y ,B z )e/mcu! L , 

where u>l and k^ are circular frequency and wave- number 
of the light wave, e and m are charge and mass of the 
electron, and c is the velocity of light. For a plane wave 
moving in vacuum in x-direction with dispersion relation 
ujl = ckL, the laser fields satisfy E y {r) — B z (t) with 
propagation coordinate r = t — x. For an electron moving 
along x(t) , we have dr/dt = 1 — f3 x . With these rela- 
tions, the equations of motion can be written with r as 
independent coordinate in the form 



■J-{l-Px) = E x, 

d Py P 



(4) 
(5) 



Integration gives the functions k(t) = 7 — p x and p y (r), 
and, recalling 7 2 = 1 + p\ + p 2 , we obtain 

1+Pl 



7(r) 
Px(t) 



2k, 
l+p| 
2k 



K 

2' 



(G) 
(7) 



Making use of k = 7 — p x = 7(1 — (3 X ) , the electron 
trajectories x(t), y(r) follow from 



dx/dr = Px/k, 
dy/dr =p y /n, 



(8) 
(9) 



in parametric form with time t(r) — r + x(t) . 

Let us consider solutions for a single cycle light wave of 
the form E y = ciq sinr, limited to the interval < t < 2ir 
and with no longitudinal field, E x = . In this case, the 
solutions of Eqs. ([4} and ([5]) and are 

Py( T ) = PyO - «o(l - cost), 
k(t) = 7(t) -p x (r) = Ko, 

where kq and p y o are integration constants. We now dis- 
cuss three important cases. 



2.2 Acceleration of electron initially at rest: kq = 1 

An electron initially at rest with p X Q = p y Q = and 70 = 1 
corresponds to kq = 1. In this case the explicit solution is 



Py( T ) = -Oo(l 
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Fig. 2. Analytic solution of longitudinal momentum p x and transverse momentum p y as function of time for an electron initially 
at rest and driven by a single-cycle laser pulse with ao = 5. Electrons at different depth in a layer are considered experiencing 
different constant longitudinal electric fields E x ; E x = corresponds to the right-hand surface of the layer, E x — 0.5 to the 
middle, and E x = 1 to the left-hand surface facing the ions of the graphene layer. For units see text. 




Px(t) 



x(t) 



-(1 - cost) 5 



7M = l + y(l-cosT) 2 , 

y(r) = a (cost' — l)dr' = — a (r — sinr), 
Jo 



2 Jo 



r' - Ifdr' 



aL3 „ . 1 . „ , 

— — ( — t - 2 suit + - sm2r), 
2 v 2 A ' 



t(j) = T + x(t). 



2.3 Collision with relativistic electron: k >> 1 



Let us now consider an electron, moving opposite to the 
direction of the light pulse and having initial momenta 
p x i = — 2a 2 ,, p y i = +2ao, and energy 71 = 1 + 2a 2 ,. This is 
equivalent to the case considered before at the time, when 
the electron has peak energy after interaction with the first 
half-wave. At this moment, it is hit head-on by the second 
pulse running in opposite direction. We call this second 
pulse the switch pulse, because it almost instantaneously 
changes p y (and thereby the phase between electron and 
wave) , but leaves p x and 7 almost unchanged. This type 
of interaction corresponds to K\ = 71 — p x \ = 1 + 4a 2 , ^> 1 
for ao > 1. Placing the electron again at the origin x\ — 
yi = at time r = 0, we obtain the solution 



For r = 7r , the electron reaches the maximum energy 
of 71 = 1 + 2dQ at time t± = tt + x(ir) = n(l + 3oq/4), 
and momenta and positions then are p x \ = 2oq, p y i = 
—2ao, x\ — (37r/4)aQ, and y\ = — nao- Momenta p x and 
p y are plotted in Fig. 2 as function of time. It is seen 
that they grow while interacting with the first half of the 
driving wave (0 < r <= tt); afterwards, in the second 
half, the particle decelerates again and comes to rest at 
t = 2tt. There is no net energy transfer between a particle 
and the light wave, unless additional fields interfere. This 
corresponds to the Lawson- Woodward theorem [5D]. For 
a clear interpretation of this theorem, one should refer to 
Esarey et al. [21] . In order to continue acceleration also in 
the second half, one has to shift the transverse momentum 
p y at r = tt from negative to positive values. In principle, 
this is possible by interaction with an additional counter- 
propagating laser pulse and will be discussed next as a 
second interesting solution of the basic equations. 



Py( T ) 
Px(t) 

7 (t) 

y( T ) 

x(t) 
t(r) 



a (l + cost), 

1 + a§(l + cost) 2 l + 4ag 



2(1 + 4a 2 ) 



1 + 



2(1 + 4a 2 ) 2 ' 

a (T + sin r)/(l+4ao), 
(l + 3a 2 ,/2)T + 2a 2 ) sinT + (oq/4) sin2T) 



2(1 + 4a 2 ) 



2\2 



x(t), 



valid for < t < 27r. Apparently, the dynamics are very 
different from the first case. Only the transverse motion 
scales oc ao as before, while all changes of longitudinal 
quantities are small of order 1 . During the first half cycle 
(0 < t < tt), transverse momentum p y turns from 2ao 
to , while energy and longitudinal momentum change 
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as Aj = Ap x = 2(Zq/(1 + 4a 2 ,) ps 1/2. This takes place 
in the short time and space interval Ax ps At ps 7r/2 
corresponding to a quarter cycle. A half-cycle laser pulse of 
same strength as the drive pulse, but running in opposite 
direction, is therefore sufficient to switch the phase of the 
electron from deceleration back to acceleration. In case 
the switch pulse extends over a full cycle, the situation of 
p y ps can be realized only for short moment. The action 
of both half- and full-cycle pulses will be compared with 
ID-PIC simulation in section 3. 



2.4 Boosting energy of a relativistic electron: |k | << 1 



A third, very instructive solution of Eqs. (Q} to ([9]) is 
obtained for \kq\ << 1. Here we start with an electron 
that has already been accelerated by the first half-wave of 
case 1 and has interacted with the half-wave switch pulse 
of case 2. Initial momenta and energy now are 



Py2 = 0, 



1 



Px2 



12 



, (l + 4o§) 

2(1 + 4a 2 ) 2 : 

1 , (l + 4ag) 



2(1 + 4a 2 ) 

K2 = 72 - Px2 = 



1 



(l + 4ag)- 



Entering with these values into the negative half-wave of 
E y = ao sin r at time r = 7r , the electron moves in the 
time interval t: < t < 2n according to 



p y (r) = a (l + cost) 
Px(t) = 



(l+4a2)(l +a 2(l +cosr )2) 



1 



7 (r) = 



2 2(1 + 4a 2 ) : 

(l + 4a 2 )(l + a 2 (l + cosT) 2 ) 1 

2 + 2(1 + 4a 2 )' 



For r = 2-7T , we find 72 ~ p X 2 ~ 8ag , and the leading 
orders of the shifts in x and y direction are Ax ps 12aQ7r 
and Ay ps 4oq7t, always supposed that ao >> 1. Choosing 
a laser pulse with ao = 5 and wavelength = 0.8 /zm, 
the single-cycle sin-pulse accelerates the electron to an 
energy of W — jmc 2 ps 2.5 GeV. For this it needs a 
distance of Ax ps 6agAL ~ 7.5 cm and a time interval of 
At ps Ax/c ps 250 ps. 

We show these results here as an interesting applica- 
tion of the one-dimensional relativistic equations of mo- 
tion. Of course, the example is difficult to realize experi- 
mentally for several reasons. It is not clear how to produce 
and transport the half-cycle switch pulse. It has a spec- 
trum with a non-zero dc component and cannot be trans- 
ported over a longer distance due to strong diffraction. An 
alternative method to switch the transverse momentum 
p y — » would be a local static B z field. This would re- 
quire \p y \ = a z k^D = 2ao, where D is the thickness of the 
magnetic layer and a z = B z /Bq with B$ = mcojL/e ps 10 4 



T'/xm/Ai. Since the switch layer should not be thicker 
than a few Xl, we would need a magnetic field in the or- 
der of B z ps £?oao/10 ps ao kT, which is quite demanding, 
but may be not impossible. 



2.5 Charge separation and longitudinal fields E x 

So far, we have considered solutions for single electrons 
only. In the case of a solid-density layer and the expulsion 
of all electrons, large longitudinal fields E x will emerge 
due to space charge and will modify the solutions given 
above. For a plane foil of thickness d and electron density 
n e , the electric field between the layers of the expelled 
electrons and the ions left behind is given by 



£ = Nk L d, 



(10) 



where Eq is the electric field in units of Eq — mcuj^/e and 

N = n e /n cnt = ujp/uj 2 L (11) 

is the electron density, normalized to the critical density 
n crit — nmc 2 / (cXl) 2 — 1.1 x 10 21 cm~ 3 /A| [fim]. Elec- 
trons, initially located at x — xq inside the electron layer, 
experience the electric field 



E x = £ (1 - xo/d), 



(12) 



which is largest for xq — on the laser-irradiated front 
side of the foil and vanishes at xq = d on the rear side. 
Neglecting the initial time period of electron expulsion 
from the foil, which is short for ultra-thin foils, the field 
E x is constant in time. This holds as long as the total 
charges to the left and to the right of xq stay the same, 
i.e. sub-layers do not change their relative order. We can 
then integrate dn/dr = E x and find 



k(t) = 1 + £ (1 - x /d)r. 



(13) 



The full solution is obtained by integrating Eqs. ([3]) - ©. 
In Fig. 2 we have also plotted the curves p x (t) including 
the effect of E x for the surface on the laser side (xq jd = 
0) and the middle interface (xo/d = 0.5) of the electron 
layer. For graphene the fields are E x — 1 and E x = 0.5, 
respectively. It is seen that the drag-back of the ions is 
very strong when compared to the leading front of the 
layer at Xo/d = 1. Only at early times the layer moves as 
a whole, but then electrons facing the ions lag behind and 
even turn around and fall back to the ion layer. 

It should be noticed that the evolution of the electron 
layers is still more complex. Almost quantitative agree- 
ment with ID-PIC simulations is obtained [23] when ap- 
plying the full analytical model of Kulagin [23] . 



3 Simulation 

3.1 Acceleration of ultra-thin electron layer 

We have performed one-dimensional particle-in-ccll (1D- 
PIC) simulations to investigate thin layer acceleration and 
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Fig. 3. Results of ID-PIC simulation: (a) Snapshot of laser field E y and electron density n e after one laser cycle when the 
single-cycle pulse has just passed the graphene layer located at x — 0. Electrons are seen as a thin layer surfing on the first 
half-cycle and marked by sharp density spikes at the surface facing the drive laser, (b) Values of 7 and j x of electrons, located 
at the right surface of the layer, are plotted as function of t. 




Fig. 4. Thomson backscattering of a semi-infinite probe pulse by the electron layer shown in Fig. 3. The weak probe pulse 
with a p = 0.001 is polarized perpendicular to the drive pulse, and its sharp front reaches the foil at t = 0. (a) Signal versus 
time as seen by an observer at 7\l behind foil, (b) corresponding spectrum, which shows a sharp cutoff at lj/ojl ~ 4yl « 36, 
well corresponding to the maximum value of j x « 3 in Fig. 3b. 



Thomson back-scattering. For illustration, we choose as 
before a single-cycle laser pulse E y = a>o sin r with < 
r < 27r, a = 5, and a layer with e = 1, correspond- 
ing to graphene. A snapshot of electron sheet acceleration 
is shown in Fig. 3a and the temporal evolution of j(t) 
and J x (t) = 7/(1 + Py) 1 ^ 2 in Fig. 3b. All electrons are 
blown out of the foil. The electron density distribution 
is plotted here three laser cycles after the laser front has 
first interacted with the foil. The thickness of the elec- 
tron layer is set to d/X^ = 0.001 initially and expands to 
Ax/Xl ~ 0.06. The density profile shows a pronounced 



peak at the surface facing the drive laser. It surfs on the 
laser wave and has substantially depleted the front part of 
the wave due to energy transfer to electrons. The electrons 
move close to velocity of light, but lag somewhat behind 
the light front at x/Xl = 3. Electron 7-values correspond- 
ing to the density spike are seen to rise to a maximum of 
7 w 19 and then fall back to unity. Though this behavior 
is qualitatively similar to that of the analytical solution 
displayed in Fig. 2, the peak value is far less than the 
analytical maximum of 7 = 2<2q = 50. The deviation is 
attributed to modifications in the radiation field due to 
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Fig. 5. Same as the case of Fig. 3, but applying in addition a full-cycle switch pulse that hits the accelerated electron layer 
approximately at t « tl- The plots show (a) momenta and (b) 7 and ~f x as function of time. The switch pulse was chosen such 
that transverse momentum p y — > for a short moment and, correspondingly, 7^ — ► 7. 



t 
E 

w 



0.4 
0.2 
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(b) 
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t/x 
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co/co 



Fig. 6. (a,b) Thomson signal and (c) spectrum of probe pulse corresponding to the case of Fig. 5 with full-cycle switch pulse. 
Notice the two conspicuous spikes in the backscattered signal, zoomed in (b). They generate the high-frequency part in the 
spectrum with power- law decay up to lo/lol ~ 500 and a shoulder extending up to lu/ujl ~ 2000. See text for origin of this 
behavior. 



energy depletion and self-radiation of the layer. Both ef- 
fects are not taken into account in the present analytical 
treatment. 



3.2 Thomson backscattering 

We now add a weak probe pulse to the simulation above. 
It propagates antiparallel to the drive pulse and is polar- 
ized orthogonal to it so that the fields of probe and drive 
pulse can be clearly separated in the simulation. The probe 
pulse is sinusoidal and semi-infinite with a sharp front, 



which first hits the foil at t = 0, simultaneously with the 
drive pulse. The back-scattered signal is shown in Fig. 
4a, as detected by an observer located at x/Xl = 7. The 
scattered pulse is strongly chirped due to the acceleration 
of the electron layer. Notice that it is also strongly com- 
pressed. It vanishes after « 8 due to density decay of 
the layer. The corresponding spectrum in Fig. 4b shows a 
broad plateau with a sharp cut-off at uj/ujl ~ 40. This cor- 
responds nicely to jx^iax ~ 3.2 (see Fig. 3b) and the rel- 



ativistic Doppler shift uj/lul 



7a;,max(l + Px,max) 



39. 



Apparently, only the velocity component [3 X in the nor- 
mal direction of the electron layer and the corresponding 
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7a; = l/y/l — (3% = j/yl + Pij is relevant for the Doppler 

shift. The full shift of uj/ujl = 47 2 lax sa 1400, expected for 
7max ~ 19 (see Fig. 3b), is suppressed by the transverse 
momentum p y of the electrons in the layer. Of course, in 
the present ID-PIC simulation accounting for three mo- 
mentum coordinates (p x , p y , p z ), but only one spatial 
coordinate x, scattered radiation is obtained only in in- 
direction. But also in 2D/3D PIC simulation, coherently 
scattered radiation will not be emitted into the oblique 
direction of electron motion. This is because of the plane 
symmetry of the coherently accelerated electrons and de- 
structive interference of radiation scattered into directions 
outside a very narrow cone in normal direction. Incoher- 
ently back-scattered radiation may be observed in other 
directions, in particular, when the density of the layer has 
decreased strongly, but this incoherent radiation scales 
with the number scatterers N sc rather than N 2 C for co- 
herent scattering and is therefore much weaker. 

3.3 Reflected attosecond pulses after single-cycle 
switch pulse 

In order to see the full 47 2 shift, one has to turn the 
transverse momentum p y — ► by some additional interac- 
tion. Such interaction can be provided by a short counter- 
propagating switch pulse. As we have discussed in section 
2.3, this pulse has to have similar strength as the drive 
pulse and same polarization. Here we present a simula- 
tion in which a single-cycle switch pulse, having ag = 2.5 
and propagating in (— ^-direction, hits the electron layer 
at time instant 4/tl = 1, when the 7-factor of the elec- 
trons is w 11.5 as seen in Fig. 5b. We choose this early 
interaction time, because it is found that the layer rapidly 
expands in this case and that the coherently scattered 
signal would be significantly weaker, had we waited until 
t/ri = 5, when 7 max ~ 19 is reached. 

In Figs. 5a and 5b, the effect of the switch pulse is 
demonstrated for electrons located at the right surface 
of the electron layer (see Fig. 3a). The transverse mo- 
mentum makes a sudden excursion to p y w 0, while the 
longitudinal momentum p x is hardly affected. This corre- 
lates to a sudden jump in "f x to almost j x ss 7 rs 11.5. 
We expect that this interaction causes a burst of short 
wavelength radiation from the coherently scattered probe 
pulse. Indeed, we observe two spikes of attosecond du- 
ration in the reflected wave (see Fig. 6a and zoom in 
Fig. 6b), a smaller one at t/rt = 7.285 and a promi- 
nent one at t/rt = 7.345. Correspondingly, the scattered 
spectrum (Fig. 6c) exhibits strongly enhanced emission 
in the frequency range oj/ujl ~ 40 — 2000. Apparently, 
the two spikes, separated by a time interval of At/r^ ss 
0.06, originate from scattering at front and rear side of 
the electron layer (see Fig. 3a). They are separated by 
a distance Ax/Al « 0.06. Actually, the two signals in- 
terfere and produce the conspicuous oscillations seen in 
the high-frequency part of the spectrum with Auj/ujl ~ 
1/0.06 » 17. This interference documents the coherence 
of two back-scattered signals, emitted from electron inter- 
faces at different depth in the layer. 



Similar to the harmonics spectra arising from reflection 
from plasma surfaces |13j . the high-frequency spectrum 
falls off according to a power law up to lo/lul as Aj 2 as 
450, consistent with j x « 7 w 11 at switch time. Also 
the spectral part beyond uj/lol ~ 450 is characteristic for 
attosecond bursts, showing no clear cut-off. The power law 
exponent is approximately —2, i.e. the fall-off is smaller 
than for the value —8/3 observed for surface harmonics 
[12113) . We attribute the smaller exponent to the fact that 
p y {t) is just touching the p y — axis in Fig. 5a rather than 
intersecting it, as it does in the case of surface harmonics 
(compare Pukhov et al. [26]). 

3.4 Relativistic mirror after half-wave switch pulse 

Let us finally discuss how a half-cycle switch pulse would 
change the results. Again we apply an additional pulse 
E = E s sin(r — r s ) with 0<T = t + x<ir, the phase 
t s chosen such that it hits the electron layer at t/rz, ~ 1, 
having same polarization as drive pulse and an amplitude 
E s taken such that p y — > is obtained approximately. 
The results depicted in Fig. 7a and 7b show that we in- 
deed find p y » and 7^ — > 7 ~ 11 for t/rt > 1 when 
selecting the appropriate parameters. Probing this rela- 
tivistic mirror as before by a perpendicularly polarized 
probe pulse, we find the back-scattered signal of Fig. 8a 
with the spectrum of Fig. 8c. On top of the reflected sig- 
nal without switch pulse (compare Fig. 4a), we observe a 
high-frequency feature (zoomed in Fig. 8b), now showing 
many fast oscillations after the switch pulse has turned 
the electron momentum into normal direction. The zoom 
in Fig. 8b exhibits two regions, low- amplitude fast oscil- 
lations for 13.285 < t/r^ < 13.345 corresponding to the 
front interface of the electron sheet and high-amplitude 
chirped oscillations corresponding to the prominent den- 
sity spike at the rear side of the layer, which is already de- 
celerating. Since the probe pulse now sees the full 7 « 11 
over an extended time, the reflected spectrum exhibits a 
clear peak at uo/lul ps 47 2 s» 480. 

4 Conclusions 

In this paper, we have explored coherent Thomson back- 
scattering from relativistic electron layers, blown out from 
ultra-thin solid foils by high-contrast few-cycle laser pulses. 
For this to happen, the laser amplitude clq has to be larger 
than the electrostatic field eo building up when separat- 
ing all electrons from ions in the foil. It requires ultra- 
thin foils: already a mono-atomic carbon layer (graphene, 
used here as reference material) produces eq ~ 1. A sec- 
ond requirement is ultra-high contrast of the laser pulses 
to avoid premature foil destruction. Though demanding, 
these conditions may be met in near- future experiments. 

Here we have investigated the generation of relativis- 
tic electron layers by schematic single-cycle laser pulses 
with ao — 5, making use of well known analytical theory 
and 1D/2D PIC simulation. The ultra-thin foils are trans- 
parent to light, and electron acceleration occurs by direct 
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Fig. 8. (a,b) Thomson signal and (c) spectrum of probe pulse corresponding to the case of Fig. 7, but with observer position 
moved to x — 13Al. The backscattered signal now exhibits fast oscillations in the time interval 13.25 < t/rt < 13.40, shown 
as zoom in (b). The small amplitude oscillations for 13.275 < i/Vi, < 13.335 correspond to the right edge of the electron layer 
(see Fig. 3a). It acts as a relativistic mirrors moving with j x ~ 7 ~ 11 and generates the sharp peak at uj/ul ~ 4~/ 2 ~ 484 in 
the spectrum. The larger amplitude oscillations, seen in Fig. 8b for 13.335 < t/rt < 13.375, correspond to the density spike at 
the left layer surface. This peak is already decelerating, and this causes the negative chirp. 



action of the laser fields, admitting an approximate analyt- 
ical description in terms of single-electrons in a plane light 
wave. Equivalent electrons then move in planes, and co- 
herent backscattering of probe radiation occurs according 
to normal mirror rules due to phase selection. Conditions 
for coherent backscattering are discussed in a companion 
paper (paper II ), and the scaling of reflectivity with foil 
density, thickness, and 7-factor are derived there [22]. 

An important point is that, due to the laws of direct 
laser acceleration, high-7 electrons necessarily have high 
transverse momentum p y and move under an angle rel- 



ative to the direction of laser propagation, which is also 
the normal direction of the electron layer. The relativis- 
tic Doppler shift of reflected light is therefore reduced to 
uj/lu l « 4 7 2 = 4 7 2 /(l +p 2 y ), since only lx = l/^/l -0* 
related to normal cc-direction is relevant. Methods to switch 
p y — > by an additional counter-propagating laser pulse 
or by a static magnetic field have been discussed. Various 
cases are demonstrated by ID-PIC simulation. 

Probe radiation backscattered from a driven electron 
layer typically shows a reflected signal with a broad chirped 
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spectrum and a cut-off at lu/ujt, ~ ^lx- The chirp is due 
to layer acceleration and may well extend over several oc- 
taves. Disrupting the smooth layer motion by a counter- 
propagating single-cycle switch pulse changes the spec- 
trum of the back-scattered probe radiation significantly. 
Having an appropriate amplitude, the switch pulse turns 
p y — > for a short moment, which is much shorter than 
a laser cycle. The probe radiation then experiences the 
full 47 2 Doppler shift for this short moment, and a single 
attosecond burst of high-frequency radiation is backscat- 
tered. A broad spectrum in the region of lu/ujt, ~ ^J 2 
is observed, which looks quite similar to the surface har- 
monics, recently observed experimentally [12j . We have 
also investigated the more speculative case of a unipolar 
half- wave switch pulse which may turn p y — ► for a longer 
time interval. The PIC simulation then shows an extended 
period of full Aj 2 Doppler shift and, correspondingly, a 
strong peak at uj /lot, Aj 2 in the spectrum. 

Though the long term goal is to use the scattered probe 
radiation for flexible generation of controlled attosecond 
VUV- and X-ray pulses, the immediate application will be 
to analyze the electron layer and its angular and temporal 
evolution in all detail. Already the few cases shown in 
this paper illustrate the abundance of detailed information 
available from the spectra. The present study is only a first 
step. Many questions arise concerning what happens with 
more realistic few-cycle pulses than those used here. Also 
mono-atomic foils may not be available immediately, and 
acceleration of foils a few nanometer thick needs to be 
studied. This is planned for subsequent investigations. 
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